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We consider a free quantum scalar field satisfying modified dispersion relations in curved 
spacetimes, within the framework of Einstein- Aether theory. Using a power counting anal- 
ysis, we study the divergences in the adiabatic expansion of (if)"^) and (T^j/), working in the 
weak field approximation. We show that for dispersion relations containing up to 2s pow- 
ers of the spatial momentum, the subtraction necessary to renormalize these two quantities 
on general backgrounds depends on s in a qualitatively different way: while {(f)'^) becomes 
convergent for a sufficiently large value of s, the number of divergent terms in the adiabatic 
expansion of (T^^) increases with s. This property was not apparent in previous results for 
spatially homogeneous backgrounds. 

PACS numbers: 04.60.Bc, 04.62.+V, ll.lO.Gh 

There are theoretical and phenomenological motivations for studying quantum fields with mod- 
ified dispersion relations (MDR) in the effective field theory framework for semiclassical gravity 
. For example, in some inflationary models the expansion of the universe can 
last sufficiently long that length scales of cosmological interest today could have been sub-Planckian 
(or smaller than a critical length for which new physics could show up) at the beginning of inflation. 
In such models, the inflationary predictions could be affected by such unknown physics, giving rise 
to so-called trans-Planckian effects 511 . A similar argument in the context of black hole physics 
has motivated the analysis of the robustness of the Hawking effect [6|] under the modification of 
the dispersion relation of the quantum fields. It has also been argued that trans-Planckian effects 
could be relevant in the astrophysical context, as for example in the physics of ultra high energy 
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cosmic rays [J]. 

In some models, the inclusion of fields satisfying MDR locally breaks Lorentz invariance. Nev- 



ertheless, one can include MDR while preserving a general 



a. 



y covariant metric theory for gravity by 



In this theory, the general covariance 



working in the framework of the Einstein- Aether theory 
is preserved by introducing a dynamical vector field called the aether field, which is constrained 
by means of a Lagrange multiplier to take a non-zero timelike value, u^n^ = — 1. 

In the semiclassical approximation, a study of the consequences of changing the dispersion 
relation should involve a proper treatment of the divergences that appear in the expectation values 
of observables associated to the quantum fields. For example, a careful evaluation of the expectation 
value of the stress tensor is necessary to evaluate the backreaction of the quantum fields on the 
background metric, since this expectation value is the source term in the semiclassical equations 
for the metric. 

The renormalization procedure for quantum fields satisfying the standard dispersion relations in 

0,Q, 



curved backgrounds is well established 



Indeed, there are well known covariant methods 



of renormalization that can be implemented in principle in any spacetime metric. When applied 
to the expectation value of the square of the field (0^), or to the mean value of the stress tensor 
{Tfj^u), one can obtain the associated renormalized quantities by making the subtractions: 



/rcn 



(la) 



(V)rcn = (T^u) - (T^.)(°^.. - (lb) 

where a superscript 21 denotes the terms of adiabatic order 21 of the corresponding expectation 
value (i.e., the terms containing 21 derivatives of the metric). For the usual dispersion relation, in 
n dimensions the subtraction involves up to the adiabatic orders 2zJ^^^ = 2int(n/2 — 1) for (cf^) 
and 2j^^2, = 2int(n/2) for the stress tensor, where int(a;) is the integer part of x (the superscript 
u stands for usual). 

In the case of MDR, the renormalization procedure have been investigated in some particular 



cases 



Id, Q, Q. For scalar fields propagating in an n-dimensional spatially flat Robertson- 



Walker (RW) spacetime, in Ref. [12t | we have considered the extension of the adiabatic subtraction 



scheme based on a WKB expansion of the field modes [15j] to the case of generalized dispersion 
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relation. There, the WKB expansion of the stress tensor was computed up to the fourth adiabatic 
order. With the use of dimensional regularization, it was shown that these adiabatic orders can be 
absorbed into a redefinition of the gravitational bare constants of the theory (i.e, the cosmological 
constant, the Newton constant and the coupling constants associated to terms quadratic in the 
curvature in the gravitational action), which correspond to the counterterms also required for the 
usual dispersion relation. However, this simple result is a peculiarity due to the symmetries of the 



spatially flat RW metric. Indeed, in Ref. 



131] we have followed the same approach for the case of 



four-dimensional Bianchi type I metrics. Restricting the calculation up to the second adiabatic 
order, we have shown that new counterterms do appear, which involve the aether field in addition 
to the metric. For instance, a term proportional to (V^n^)^ in the aether Lagrangian is needed in 
order to absorb the divergences in (T^^)'^'^'^ (in addition to the usual Einstein-Hilbert action). The 
point is that in a spatially flat RW background these new counterterms are indistinguishable from 
the usual ones. Concretely, once evaluated in this background, the stress tensor obtained from the 
variation of the most general action for the aether field containing two derivatives, turns out to be 
proportional to the Einstein tensor. As there are strong constraints on the parameters associated 
to terms containing two derivatives of the vector field [al, the new counterterms of second adiabatic 

n 

order should be carefully chosen to make the theory consistent with observation |13l |. 

In both spatially flat Robertson- Walker and Bianchi type I backgrounds, given a MDR such 
that the frequency behaves as a; ~ lA;!'* for large values of \k\, one can show that divergences appear 



up to 



12, 



13i 



2zL, = 2int(^--j, 2jL, = 2intf- + ^), (2) 



where the superscript h stands for (spatially) homogeneous. Note that the correct values for the 
usual dispersion relation are recovered for s = 1. Eq.Q suggests that the higher the power s of \k\ 
the milder would be the divergences. For example, in four dimensions, this result indicates that if 
s > 4 the divergences of the stress tensor are contained only in the zeroth adiabatic order. However, 
as will be shown below, this is valid only for the spatially homogeneous backgrounds considered so 
far. 

In this letter we will consider weak (but otherwise general) background fields, and use a power 
counting analysis to investigate up to which adiabatic order the subtractions in Eq.([T|) have to be 
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carried out in order to get finite, physically meaningful expectation values. In the case of interacting 
quantum fields in Minkowski spacetime, higher spatial derivatives improve the UV behavior of 
Feynman diagrams (see for instance Here, we will show that while such improvement also 

occurs for ((/>^), this is not the case for {Tf^u)- 

By adopting the semiclassical approximation we consider a massive quantum scalar field (j) 
propagating in a classical curved spacetime and coupled to a classical aether. The classical action 



for the scalar field is given by 



17 



(3) 



where n is the spacetime dimension, g = det{g^^), is the standard Lagrangian of a free scalar 
field 

1 



(4) 



with R the Ricci scalar, Ccor is the corrective lagrangian that gives rise to a generalized dispersion 
relation 



<i,p<q 



(5) 



with =-L^ Va(-L7 V^c/)), where -L^i/= g^u +u^Ui, and is the derivative operator associated 
with g^y. 

After varying the action ^ with respect to one can write the equation for the associated 
Green's function G(x, x') as 



□ - (m^ + iR) - bgp {H'^V'^P + RPV^'^) 



G{x, x) 



--5{x — x), 



(6) 



where H is an operator given by 



(7) 



with a° = u^'Vf,u"'. 



In the weak field approximation, we consider the linearized Eq.([6]) about the flat background 
solution with Minkowski metric rj^i, and constant aether [l^. Adopting Minkowski coordinates 
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{x^, ^) such that = 5^, the fields can be expanded as 



Here and in what follows, we keep only linear terms in h^i, and w^j. The constraint u^u^ 
yields = h?^ /2. After choosing the Lorentz gauge 



5^ V = 0. 



hinj hi,. 



(8a) 
(8b) 

= -1 

(9a) 
(9b) 



with h = hp, we obtain 



g\a = {rj>'- - e'')d^d, 



(10a) 



<i,p<q g,p<q q,p<q 

where A is the Laplacian and /C is an operator of first order in the perturbation fields h^i, and 
v^. The explicit expression for the operator /C can be easily found in the particular case in which 
the perturbation fields depend only on the time coordinate x^. In the general case, although more 
complex, it can also be worked out. However, the expression in Eq. ()10bp will be enough for our 
present purposes. 

In order to solve Eq.(l6|), we first split G = + G^, where the superscripts refer to the order in 
/i^j/ and v^. The zeroth order propagator G^ satisfies 



-^0^ + A - - 2 ^ bgpA'^+P 
q,p<g 



Hence, the corresponding Feynman propagator is 



G%{x,x') 



G^x,x') = -6{x-x'). 



(11) 



(2vr)" [-ko^ + u;^\k\^) - ie] 



(12) 



where kx = ij^^k^x'^, and 



2n 1,12 



q,p<q 



(13) 
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with k = {ki, kn-i). 

Then, the first order correction to the propagator satisfies 



-a^ + A - - 2 ^ b^pA^+P 



and the corresponding solution is 



G^{x,x') 



h^'d.d^ - + ^(u;2(-A) + A) + /C 



(14) 



G\p{x, X ) 



d^yG^x,y) 



h^''{y)d,d, - )-t%d,d,h{y)) + -^{u;'{-A) + A)) + /C(y) 



G%{y,x') 
(15) 

In what follows, for the sake of convenience, we work in Euclidean spacetime with kn = —ik^ and 
. The Euclidean Green's function G^; is given by Gp = iGe- 
For our power counting analysis we consider, as an example, only one contribution to Gg, which 
we choose to be 



gU^,x') = - I ryhiy)G%{x,y)Lo\-A)G%{y,x'). 



(16) 



One can readily check that this term, along with the one proportional to hP^{y)dodo and similar 
terms within JC{y), are the most divergent ones. By replacing the zeroth order solution G^ into 
this equation, introducing the Fourier transform of h, 



h{p) = / d^ye-'Pyh{y), 



(17) 



and performing some trivial integrations, we find 



„iA:(x-x')^2(|^|2 



(18) 



In order to derive the superficial degree of divergence of each adiabatic order we start by ex- 
panding the integrand in Eq. (jl8p in powers of Pi {i = l...n), 

.K^ _ 1 f ^> Jpxuf^\ f d^^k e**=(^-^')a;2(|fc|2) f 1 



qpix,x) = - I - — ~—e , , , _ , 

2j (2^)- "''J (27r)- [fc2+^2(|^|2)]2 Vl + ep 



1 f d> , 



2 J (27r)' 



(27r)- [fc2 +^2(|^|2)]i 



+ 00 



(19a) 
(19b) 



r=0 
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where 



2knPn+pl+C0^{\k+p\'^)-UJ^ 



(20) 



The expectation value 
can be decomposed as 



ki+u^im 

b"^) is given by the coincidence Hmit of ReG_E(= ImG_p). In such hmit, 





J (27r)" 


Ap) = J 


d'^-'^k 


(27r)"-i 


I{pAk\)= j 


dkji \ 

r=0 



h{p)T{p), 



2)]2 



(21a) 
(21b) 

(21c) 

do not con- 



Notice that, due to the remnant symmetries, the terms with odd powers of Pn or 
tribute, so there will be no odd adiabatic order contributions. 

Prom Eg. pi I) it is possible to derive the superficial degree of divergence of X(p) for each term 
with a given power of pi. However, it is more instructive to compute explicitly the integral in 
Moreover, in order to appreciate better the difference between the cases in which the background 
fields depend or not on x, we analyze them separately. 

For background fields that depend only on x„, we replace in Eq. (j2ip the n-dimensional Fourier 
transform h[p) by h[p) = {2tt)^~^ 6^~^ {p)h{pn) . Performing the integral in A;„ of Eq. (j21cp it is easy 
to see that 

2i 



iiP,\k\) = Y.' 




(22) 



.=0 ^'m' 

where ai is independent on k. Therefore, substituting Eq. (j22p into Eq. ()21bp . by means of power 
counting, we find that when the MDR is such that u;(|A:p) ~ for large values of \k\, convergence 
of the 2z— adiabatic order of I{p) is guaranteed if —s{2i + 3) + 2s + n — 2 < —1. That is, given the 
highest power s of in the MDR and the spacetime dimension n, we expect that the maximum 
adiabatic order of {(f)'^) that contains divergences will be given by 

1 



2/' 



2int 



n 



2s 



(23) 
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This result coincides with the one given in Eq.([2]). It is worth to remark that the same result can 
be obtained using the weighted power counting analysis described in Ref. [l^ . 

For background fields that depend only on x, we introduce the spatial Fourier transform h{p) 
related to h{p) by h{p) = {2TT)6{pn)h{p). In this case, performing the integration in kn we obtain 



where Pr does not depend on k. This equation should be compared with Eq. (j22p . In this case, the 
odd powers in \p\ vanish once the angular integrations in Eq. (l21b|) are carried out. For large values 
of k, u;^(|fc + — u;^(|A;p) ~ • k), then power counting indicates that the 2i— adiabatic 

order will be finite if 2f > n — 1 — s. Hence, in this case, we expect that divergences will appear up 
to the 2zmaa;— adiabatic order of ((/>^), where 

2C,, = 2int( ^~^^"" ). (25) 

Here the superscript g stands for general, since it can be easily shown that this is also the superficial 
degree of divergence for arbitrary backgrounds (i.e. h/^n, = h^iy{x,Xn)), = v^{x,Xn))- Note that 
while the values of 2imax given in this equation coincide with those obtained from Eq. (j23p for n < 4, 
it is generally larger for higher number of dimensions. For example, for n = 5 and s = 2 we have 
2i^ax = 2 while 2i^„^ = 0. 

We have performed our power counting analysis for the particular contribution to given in 
Eq. (jl6p , and one may wonder if other contributions to the propagator could cancel the divergences 
in g^. In order to check explicitly that this is not the case, let us consider a MDR of the form 
(j^d/cp) = TV? + |A;p + 26ii|A;|^ and compute the operator K, defined in Eq. ljlOp . For background 
fields that depend only on x, after performing some integrations in Eq. ijlSp . we arrive at 

1 f (F'-^p ^^.^ f d^-^k { ho 



ReG^(x, x) 



'-00 



2 J (27r)"-i J (2^)"-i (^(|fc|2)+^(|fc+^2) 

^ 2cu(|fc|2)c^(|fc+pl2)[^(|^|2)+^(|^ + ^2)] J' 

where, up to second order of pi, 
Up) = 26ii I h{p)\k\'' - hooip) [\p\^\k\^ - 2{p- kf] - 2 J2 h'^hkj [\p\^ + 2|fc|2 + 2p- k] i . (27) 
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Then, we expand the terms between brackets in Eg . (1261) up to second order in pi to obtain an 
integral expression for By using dimensional regularization, we perform integrations by 

parts and discard surface terms to express all the integrals appearing in {(j)^)^'^^ in terms of only 
two of them. In this way, we obtain: 



2> (2) _ "^-1 



8(27r)"-i i V 6 



1 3n - 7 



(28) 



where the factor Jln-i = 27r^" ^^^"^ /T[(n — l)/2] comes from the angular integration, = — A/i/2, 
RIq = -A/ioo/2, and 

("-3) 

dx — ^rr-r^ (29a) 

'^ = 1 '^im^- ^'''^ 

Here we see that while the integral /a converges in n = 5 dimensions, there appears a new integral 
/3 which is proportional to h\\ and diverges as n ^ 5. Moreover, the divergence proportional to 
i?oo is '^ot purely geometric. Indeed, it can be covariantly written in terms of the metric and the 
aether field since, to linear order in the perturbation fields, i?QQ = R^yU^vy = V^a'^. On the other 
hand, when ((^^)^^^ is computed for background fields that depend only on the time coordinate, 
the second term in Eq. (j28p does not appear and, in agreement with Eq. ()23p . this adiabatic order is 
finite. However, as cancellations are not to be expected for general background fields, divergences 
will generally appear up to the 2imaa;— adiabatic order given in Eq. (f25ll . 

The differences between the two cases we are considering are better evidenced when the power 
counting analysis is applied to the stress tensor of the scalar field. The expectation value of the 
Euclidean stress tensor of the scalar field can be expressed as the coincidence limit of a nonlocal 
derivative operator applied to the two-point function G\;{x,x'). For example, the coincidence limit 
of the derivative with respect to Xn and x'„ of G\;{x,x') will contribute to the stress tensor. For 
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our analysis let us consider the following contribution: 

^:e'P-h{p)T{p), (30a) 
(|^|2)r(p,|fc|), (30b) 





r dj'p 

J (27r)" 


T{p) = j 




(27r)"-i 


npAk\) = j 


r-=0 



(30c) 



+ a;2(|A:|2)]2 

When the background fields depend only on Xn we find that T(p, \ k\) can be written as 

+00 / \ 2i 

T{p,\k\) = y^^\^^] , (31) 

with 7j k-independent. After substituting this expression into Eq. (|30b|) we obtain that divergences 
are contained in terms with up to powers of p„, where 



2^;^,, = 2int(^- + ^J, (32) 

in agreement with Eq.(l2]). 

On the other hand, for background fields that depend only on x we obtain 

-'^^-|^(^-')' 

where C^r does not depend on k. Recalling that odd powers in |^ do not contribute and using again 
that u;^(|fe+p|^) — a;^(|A;p) ~ |A;p*~^(p • k) for large values of A;, we arrive at 

2jl,, = 2int( ^~2 + " ), (34) 

whence we see that, contrary to the previous case, Ijmax increases with s. For instance, for n = 4 
and s = 1, 2 divergences appear up to the fourth adiabatic order, while for s = 3 the sixth adiabatic 
order is also divergent. Hence, we expect that for s > 3 the renormalization of the expectation 
value of the stress tensor will require counterterms of adiabatic order higher than four. 

Analogously to the case of (0^), for a general background we do not expect that divergence 
cancellations occur. Therefore, we conclude that the subtraction in Eq.([T|) should be performed up 
to the adiabatic orders 2imaa; and 2jmax given in Eqs. (f25]) and (fM|) . respectively. In particular, in 
order to renormalize the semiclassical Einstein- Aether equations it will be necessary to introduce 
all possible counterterms constructed with g^i, and up to the 2jmax— adiabatic order. 
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